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Abstract 

n 

We consider a chain of Josephson-junction rhombi (proposed originally in |1]) in quantum regime. 
In a regular chain with no disorder in the maximally frustrated case when magnetic flux through each 
rhombi Or is equal to one half of superconductive flux quantum Oo, Josephson current is due to correlated 
transport of pairs of Cooper pairs, i.e. charge is quantized in units of 4e. Sufficiently strong deviation 
5<j> = jo,, — <I>o/2| > h^^^ from the maximally frustrated point brings the system back to usual 2e-quantized 
supercurrent. For a regular chain SO'^ was calculated in \%. Here we present detailed analysis of the effect 
of quenched disorder (random stray charges and random fluxes piercing rhombi) on the pairing effect. 



I. INTRODUCTION 



Pairingof Cooper pairs in frustrated Josephson junction arrays was theoretically proposed re- 
cently UO Isl 0] in the search of topologically protected nontrivial quantum liquid states. The 
simplest system where such a phenomenon could be observed was proposed by Dou9ot and Vi- 
dal in iQ]. It consists of a chain of rhombi (each of them being small ring of 4 superconductive 
islands connected by 4 Josephson junctions) placed into transverse magnetic field, cf. Fig. 1. 
It was shown in [JJ that in the fully frustrated case (i.e. magnetic flux through each rhombus 

= j(^o = ^) usual tunnelling of Cooper pairs along the chain is blocked due to destructive 
interference of tunneling going through two paths within the same rhombus, while correlated 2- 
Cooper-pair transport survives. Sufficiently strong deviation 54> = |^»^ — ^»o/2| > 54>'^' from the 
maximally frustrated point brings the system back to usual 2e-quantized supercurrent. 

In ref. |5] rhombi chain was studied for the experimentally relevant situation when its Coulomb 
energy is determined by capacitance of junctions and not by capacitance of superconductive is- 
lands themselves. Expression for the critical deflection 64>^ was derived. However, ref. [5] dealt 
with regular chain with no disorder. In any real system two intrinsic sources of disorder are always 
present: a) some weak randomness of fluxes penetrating different rhombi (due to unavoidable 
differences in their areas), and b) random stray charges which produce, due to Aharonov-Casher 
effect, some random phase factors to the phase slip tunnelling amplitudes leading to suppression 
of quantum fluctuations in the chain (compare analogous discussion in |^). 

In this paper we adapt method used in ref. jfl for the case of rhombi chain with quenched 
disorder and study influence of random stray charges and random fluxes in rhombi on the crossover 
point = 54>'^' between 4e- and 2e-regimes. It turns out that in a long chain with large Ej/Ec the 
pairing effect is rather stable under influence of disorder. 

The rest of the paper is organized as follows: in Sec. II we derive effective Schrodinger equation 
describing rhombi chain with quenched random stray charges or random fluxes piercing rhombi; in 
Sec. Ill we discuss influence of quenched stray charges on pairing effect, derive "phase diagram" 
of the chain with fixed realization of disorder, calculate probability 'P4e{Ej/Ec,N,d'^) to find the 
chain in the regime with dominating 4e-supercurrent and finally estimate critical deflection from 
the maximally frustrated point 54>'^' destroying 4e-supercurrent in a chain with stray charges; in 
Sec. IV we consider modulation of the supercurrent in a clean and disordered chain by external 
capacitively coupled gates; in Sec. V. we analyse influence of randomness of fluxes in rhombi on 
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Figure 1 : The chain of rhombi with random stray charges. Also shown is an external capacitively coupled 
gate (c.f. Sec. IV). 

pairing of Cooper pairs. Finally, in Sec. VI. we present our conclusions and suggestions. 

II. EFFECTIVE HAMILTONIAN IN PRESENCE OF DISORDER 

We study a chain of N rhombi shown in Fig. [U Each rhombi consists of four superconductive 
islands connected by tunnel junctions with Josephson coupling energy Ej = hl^ /2e; charging 
energy Ec is determined by capacitance C of junctions, Ec = e^/2C (we neglect self-capacitances 
of islands which are assumed to be much smaller than C). Below we consider Josephson current 
along the chain of A'^ ^ 1 rhombi and assume that the chain is of the ring shape, with total magnetic 
flux 4>c inside the ring. We also denote by 4>" the flux through n-th rhombus and define phases y 
and cp„: 



A regular chain with no disorder (no random stray charges and all (p„ = cp) is described by the 
imaginary-time action 
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and additional conditions 
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f;(-eF'-ei-")=T. (4) 

n=l 



Here the variable G,!'"^ is the phase difference across the m-th junction in the n-th rhombus (see 
Fig.|D. 

This regular model was analyzed in details in ref. |5] within the limit Ej ^ Ec- MEc = the 
phases oi™^ are classical variables which do not fluctuate. Emerging classical states of the chain 
|m, {ofj}) can be characterized by a set of binary variables {of,} (one for each rhombi) and an 
integer-valued variable m. Note that each binary variable of, can be considered as z-projection of 
spin i ascribed to each rhombi. Energies of these classical states are 
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(y- TIN/ 2 - TlS^ - 2%mf - VlhS-Ej + Const. (5) 



sl = -ol, 5^=-£a;, Y = Y+^. (6) 

Finite Ej /Ec gives rise to quantum phase slips (QPS) in the chain, which mix different classical 
states leading to formation of truly quantum ground state. Let us denote as D the amplitude of a 
QPS in one contact. At large N ^ I this amplitude does not differ from the "spin flip" amplitude 
for a single rhombus at ^ 4>o/2. Li this approximation we can use result from Ref. ^\ 

V ^ k {e]Ec) '/'exp (^-1.61 . (7) 

Here A; is a numerical factor of order one which slightly depends on Ej/Ec (see ref. |Q] for details). 

In Ifl it was shown that calculation of the persistent current in the large-A'^ limit can be reduced 
to the solution of a Schrodinger equation for a particle having a large spin S = jN moving in a 
periodic cos-like potential, with appropriate boundary condition. 

More precisely, the quantum ground-state energy E of the chain in the limit Ej ^ Ec can be 
obtained from solution of Schrodinger equation 

-^ + {E- Iwco^lx- Y^Sl + lhY, = 0, (8) 

rt=l n=\ 

where 

16NE 6ANX) , 8A^6 

E = ^ , w = , h = —y- . (9) 
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Here V — amplitude for quantum phase slip in one rhombi and 5 = TT — (p. Function 
\|/ = {On})- Sf^ and Sf^ are standard operators of x and z projections of spin j acting on o„. The 
magnetic flux inside the whole ring enters the problem through the twisted boundary condition 
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^inS^+inN/l^ (x + 71/2, o) = e'Mx, o) , j = j+^ (10) 
on the basis of this formulation of the problem critical deflection S^)*^ was derived 
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(^^"K^O.li^-j 4>o (11) 

The aim of this section is to derive effective Hamiltonian similar to (EJl for the chain with 
random fluxes 4>" through rhombi or random stray charges. 

To account for random stray charges we present the electrostatic charging energy in the form 

N 



"0= t t ^[c-']:;;;(eit*-«;^'')(G&'-4?) m) 

III ,112=1 ki ,k2=l 



where \ ^^^^ is the matrix of inverse capacitances. Indices n and k numerates superconducting 

(k) (k) 

islands (3 rows, islands in a row). Qn ' — charge of the n-th island in k-ih. row. Parameters ^„ ' 
are determined by the random stray charges. 

Starting from this charging energy we may derive an additional term in action Q emerging 
from presence of random stray charges (compare with [6]) 

55 --i dtl^l^ pf-^ (13) 

(k) 

Parameters p'" can be expressed in terms of charges qn . Corresponding expressions are a bit 
cumbersome and we do not present them here. Below we write down some special combinations 
of p"^ which we will need in our paper. 

Additional term ( fT?l has a form of total derivative. Hence it does not change neither the clas- 
sical states of the chain nor the classical equations of motion and the real part of classical action 
on a single tunneling trajectory. The only effect of this term to give a tunneling amplitude along 
each path its on phase factor. Since there are several QPS trajectories between two classical states, 
all having the same real part of tunneling action, these additional phase factors may give rise to 
destructive interference of tunneling processes, leading to reduction of total matrix element, con- 
necting two classical states. 
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Following ref. \y] and taking into account complications due to random stray charges we write 
tight-binding Hamiltonian 



- V ^ ( in 1 

//|m,o>=£™oh,a>+- J^exp<^ — {3pl~ pI- pI~ ft) ^S+|m,a> + 

n=l ^ ' 

^ E {~f {^Pl-Pi -Pl-Pt)]^n\fn-I,(s> + 

f^oxp!^--^{3pt-pi-pl-pl)^d+\m-l,o>+h.c. (14) 



2 , 



Performing Fourier transformation over variable m according to 



,a) = J^exp|2/|^2m-^+5^ + ^^x| |m,a) , (15) 



we obtain the effective Schrodinger equation 

N ^ N 

E anS'l — 2wsin2;c ^ 

n=l n=l 



1^ + I £ - 2wcos2;c - 2wsin2;c bnSl + 2hS^ ) V = 0, (16) 



and twisted boundary condition 

g.7r5HmA?/2^ (;c + 71/2, o) = o) . (17) 

Here parameters E, w and are described by equation Un and bn are random coefficients which 
can be expressed in term of random charges qi^ 

— 5 t>n — (io) 



q:^pI-pI= -^tt qi^pI-p\= - 4 1 1 (19) 

Here we measure charges q'n' in units 2e. 

We turn now to generalization of (EJ) for the case of random fluxes in rhombi. Imaginary-time 
action for this problem is given by equation Q but additional conditions ^ should be changed to 

4 

£0l'")=(p„, R=l,2,..,A^, (20) 
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Following the same steps as before we see that to account for disorder in flux piercing the rhombi 
one should just replace the last term in (HI) by 

if^KSl (21) 

n=l 

Here = 8A^(7r - (p„) /n^. 

Now we have effective Hamiltonians of the chain in presence of disorder. In the rest of the paper 
we will analyse these Hamiltonians in order to find out influence of disorder on the crossover point 
between 4e- and 2e-regimes. 

ni. INFLUENCE OF RANDOM CHARGES ON CROSSOVER POINT 

In this section we study influence of random stray charges on the on pairing effect in rhombi 
chain. It is important to note that, generally speaking, even at maximally frustrated point /z = (in 
contrast to regular chain) symmetry properties of Schrodinger equation (fT6t do not prohibit 2e- 

(2) (3) 

supercurrent. This is due to the fact that asymmetric realizations of random charges with qh j^qk 
(and thus bn 7^ 0) break symmetry between two tunneling trajectories of a Cooper pair within the 
same rhombus. Of course, random charges preserve classical states of the rhombi chain and in 
a chain with no QPS there would be no 2e- supercurrent at maximally frustrated point, but full 
quantum Hamiltonian of the chain does not possess corresponding symmetry. Nevertheless, we 
will see later that for typical realizations of random charges and at = 4>o/2 the 2e-supercurrent 
is small as compared to 4e-supercurrent. The reason for that is as follows: if at least in one rhombi 
bn = { or Un = 0) then 2e- supercurrent is prohibited. Here we start with analysis of general 
situation of rhombi chain with random charges and h and then make some conclusion on 
disordered rhombi chain in maximally frustrated point. 

We investigate the grand partition function for the system described by equations dT^fTTl) 




(22) 



with P being the inverse temperature, (3 
functional of x{x) defined as 
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. Operators f/„ act in the space of spin j each and are 




(23) 
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fn (x) = a„ COS 2x(x) , gn (x) = bn sin 2x{x) (24) 

As was shown in jsl, in a regular chain with Ej ^ Ec the borderline between 4e- and 2e- 
supercurrents is at rather large — 4>o/2, in the sense that at the crossover point h ^ w. So 
in this paper we will also consider this limit only. 

Under such condition equation (l23t can be solved for arbitrary functions /(x) and ^(x). For 
Trf/„(|3) we then find 

TrC^(|3) = exp + ^ (MO) - igniO)) (/„(|3) + /^„(|3)) + 

^ / rfXirfX2(/„(Xi)/„(X2)+^.(Xi)g„(X2))e-''l^'-^^l + 

/^^^JxiJx2/„(xi)^,(x2)e-''l^'-^2lsign(xi-X2)^ (25) 

From equation (l25l) we derive effective action for variable x 

Z= /^©x(x)e-^W^)l S = Sbound + S, (26) 
Jo 
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Shound = -T72 ianCOs2x{^) + ibnSm2x{^)) {ancos2x{0) ~ibnsm2x{0)) (27) 

n=l 

Sx= dx A — / (ixi(iX2Cos2x(xi)cos2x(x2)exp(|xi — X2I) — 

Jo 2 8 Jo 

5— / JxiJx2sin2A:(xi)sin2;c(x2)exp(|xi — X2I) — 
8 Jo 

iC— JxiJx2COs2A:(xi)sin2;c(x2)exp(|xi — X2|)sign(xi — X2) (28) 
4h Jo 

Where 

A=£a2, B=f^bl C=f^a„bn (29) 

«=1 «=1 ;j=l 

Note that similar approximation was used previously in |5] for a regular chain. In a regular 
chain semiclassical analysis which do not rely on linearization with respect to w /his also possible. 
It turns out that exact value of the critical deflection d'^^ differs only by 12% from the one obtained 
by linearization even for the case when at the crossover point w/h= 1. 

Action (l28t is nonlocal and looks a bit terrific, but since 

/a' 



~ ^'^) = -2/?6(xi - X2) 



(30) 



we can perform Hubbard-Stratonovich transformation and derive a representation for grand parti- 
tion function with local action (after redefinition of time scale according to x — > t/h). 



/ 2)Jc(x)2)j(x)!Dz(x)e"^W")] (31) 
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S = h dx{ 1 \-aidycos2x + i^idysin2x+ 

Jo \ 2 2 
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a2dzsm2x + i^2dzcos2x — — (P^sin 2x + [32Cos 2x) ] (32) 



Here tti, aj. Pi and P2 should satisfy equations 



a? + p2 = 9, P? + «i = ^, aipi-a2p2 = ^ (33) 



and 



As we see from equations (l33l) we have some freedom in definitions of these parameters. 
Namely, if we set 

l~A [a 
0Ci = Y-sinKi, p2 = W-cosKi (35) 

«2 = Y^sinK2, pi = W— COSK2 (36) 



sin(Ki - K2) = ^= (37) 



then (l33t reduces to 



and we have only one equation for two parameters Ki and K2. This freedom will not be important 
for us. In present paper for the reasons described below we will mostly consider action (1^^ under 
condition C = and chose 

ai = yf, = Pi=0,p2 = (38) 

The main feature of the result presented above is that we have reduced original problem includ- 
ing large number of random variables (~ A^) to a problem with only three random parameters A, 
B, and C. More over, since A, B and C arise as sums of large number A^^ of independent random 
variables (see equation (l29b '). it is natural to expect that they obey Gaussian statistics. Let us as- 
sume that parameters and 2,^ in (fT9t are uniformly distributed on [—1,1]. It corresponds to 



strong fluctuations of stray charges from sample to sample. Under such an assumption one can 
easily find 

(A) = {B) = j 

\ / \ / 16 64 ' ^ ^ 16 64 ^ ^ 

(c) = o. = ^ 

Our strategy in the rest part of this section will be to analyse properties of the system described 
by equations (1^ [TTl) with fixed A, B and C and then make some statistical analysis. Under 
fixed A, B and C the main subject of our investigation will be whether the system is in regime of 
dominating 4e-supercurrent or not. 

First of all let us analyze action (l32b for a trajectory where x{x), y{x) and z(x) are constant. We 
have 

= ^hy-^ Vaxdycoslx + a2dz%m2x-—{^\?,\Y?2x + ^lco^^2x) \ (40) 

This static action has two groups of minima (we call them even and odd, suppose A> B) 

x^nn, y = -aid, z = (41) 

and 

x^-+Kn, y = aid, z = (42) 

where n is an arbitrary integer. All these minima correspond to the same value of Sst- So we have 
to consider two types of tunnelling trajectories. Trajectories of the first type connect minima of 
the same group, i.e. "even-even" and "odd-odd", and corresponding variation of the variable x 
between minima is ±71, whereas y returns to its original value. Trajectories of the second type 
connect minima of opposite parity (i.e. opposite signs of y), and change x variable by ±|. It 
is not difficult to see from Eqs. (ll5ll6ll7n . that increment Ax: of the variable x along tunnelling 
trajectory is in one-to-one correspondence to the elementary charge transported along the rhombi 
chain: qo = ^Ax. Therefore trajectories of the first type lead to 4e - supercurrent, whereas trajecto- 
ries of the second type produces usual 2e-quantized supercurrent. In semiclassical approximation 
amplitudes of the supercurrent components are determined primarily by the classical actions on 
corresponding trajectories: 

/(Y) = he sinY+ ke sin(2y) , (43) 
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where 

he=A4eeM-SE): he = A2eexp{-Sjf) , (44) 

To find out whether 4e-supercurrent dominates in the system, we should compare classical 
actions on the trajectories of two types described above. So we examine classical equations of 
motion for action (l32b . For simplicity here we put coefficient C to zero. It can be shown that 
relatively small C of order (compare to (l39t ) is not important for our future purpose. 

i:' + 2aiJjsin2.x: — 2a2<izcos2A: = (45) 



y — y = tti J cos 2x (46) 

z — z = CHid sin 2x (47) 

Here we will analyse equations (l45l l46l l47l) and find analytic expression for the borderline 
between 2e- and 4e-regimes under conditions 

A(f/N->\, A-5<A (48) 

We also present results of numerical computation of the borderline in general situation. 

Let us start with 2e-trajectory. Note that for variables x and y characteristic frequency is 1. Let 
us suppose that on 2e-trajectory x varies slowly so that characteristic frequency cO;,; <^ 1 . Than we 
can eliminate y and z from equations of motion in adiabatic approximation and obtain 

x-dl^mAx = Q (49) 



^i-A^ + 2j2(A + 5) ^^^^ 
We see that under conditions (l48t (Ox ~ is indeed small. Equation (l49t has solution correspond- 
ing to 2e-trajectory. 

x(x) = ^arccos(-tanh(2Jix)) , (51) 
Classical action on this solution is 



(52) 
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Let us now turn to examination of 4e-trajectories. Here we suppose that x is a fast variable i.e. 
(Ox ^ 1- We can neglect then variation of y and z on classical trajectory and put y = —aid, z = 0. 
This is consistent with boundary conditions for 4e-trajectory. Equation for x becomes 

X J^sin2x = (53) 

N ^ ^ 

Again, due to (l48t . we see that x indeed varies fast on 4e-trajectory since (Ox ~ ^/Ad^/N ^ 1 . So 
we find 4e-trajectory and corresponding action 



x{x) = arccos ( -tanh ( 2y ) ) , (54) 



S4e=4h^— (55) 

Comparing equations (l52t and S55h we find the line of crossover between 4e- and 2e- regimes: 
the set of points (A, 5) such that = S2e- 



Bd^ IfAdW Ad^ 



N \l \ Nj ' N 
The result presented above was derived under assumption C = but proceeding in the same 
way with C < one can show that at large nonzero C does not affect the crossover line (15^ . 
Therefore in the rest part of this paper we will completely ignore coefficient C. 

The borderline obtained from numerical solution of classical equations of motion (l45l l46l l47b 
and its asymptotic form (l56b is presented on figure |2| Instead of A and B we have used here more 
convenient variables 

.(A -5) d^{A + B) 

Note that by definition v > \u\/2. From (l56t we see that in terms of u and v at large v the borderline 
between 2e- and 4e-regimes is given by 

" = /(v)~^^4, v^oo (58) 
V — 2 

On figureEh "phase diagram" of the disordered chain is presented in a more intuitive manner: we 
emphasize here that 4e-supercurrent exist in a small vicinity of the maximally frustrated point. 

Equipped with this result, for any given set of quenched random charges (characterized by 
coefficients A and B) we can (in principle) determine whether 4e-supercurrent dominates in the 
chain. Of course experimentally we have no access to such quantities as A and B. Hence we need 
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■ boundary u=|v|/2 

■ L=4v/(v-2] 



SO D 
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Figure 2: "Phase diagram" of disordered chain. SoUd Unes on the both figures mark the crossover le- and 
4e-regimes. Note that these lines does not correspond to any phase transition. The crossover however is 
sharp at large N since actions SA,e and S2e are proportional to the number of rhombi. Subplot a) presents 
crossover line in variables u and v defined by dSTl which are useful for calculations. Subplot b) depicts the 
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same in a bit more intuitive way. Here we introduce factors A = ^ ^^^i j and B = (^-^rfgij describing 
realization of disorder; (SO'^),,^^ — critical deviation for a clean chain with same Ej/Ec- 

some statistical description of rhombi chain with random charges. Such a description is provided 
by the probability (P {Ej /Ec,N,8) to find dominating 4e-supercurrent in the system. 

Assuming Gaussian statistics for u and v and taking into account (l39l) one can derive probability 
distribution for u and v 



SN 



exp 



2Nu' 



exp 



32A^(v-vo) 

J4 



(59) 



vo = ^ 



Required probability can be written as 

neiEj/Ec,N,d)^l-2 dv duP{u,v) 

Jo Jo 



(60) 



(61) 



Maximum of probability distribution (15^ lies at m = 0, v = This point for J = 9 is marked 
on figureElwith a cross. At sufficiently large d we can replace /(v) in (FTT ) by 4. We then get 



2 IT 

{Ej/Ec,N,d) = 1-^ ["^^ duexp(-u^) = 1-Erf 

Jo 



I 32N 
'd^ 



(62) 
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Figure 3: Probability as function of deviation from maximally frustrated point. Subplots a), b), c), d), 
correspond to = 10, 20, 30, 70. 

Let us introduce parameter K and define critical deviation from the maximally frustrated point 
54>c as deviation under which the probability (l62b equals K. Reasonable choice for K is 0.5 or 0.75 
or something else. From (l62b we then get central result of this section 

54>.._ (Erf-i(l-K))^/^ 1 

4>o 23/2;tl/3 \Ej) ^ ^ 

Comparing (l63t and (fTTl) we see that for any experimentally reasonable critical deviation from 
the maximally frustrated point does not differ from the one in a regular chain. However one 
should remember that theory presented above has several limitations. To clarify this question 
we present here several graphics for (P/^ei^.Ej /Ec,N) at different Ej /Ec and together with 
the boundaries of the validity region of the approximation used (fig. |3l). Subplots a), b), c), d) 
correspond to different number of rhombi = 10, 20, 30, 70 respectively. Solid and dashed lines 
1,2,3,4 correspond to different Ej/Ec = 10, 12, 15,20. Dashed lines were obtained with the aid 
of equation (l62t relying on the condition > I whereas solid lines were produced by exploiting 
equation (l6TT) with /(v) determined from numerical calculations. The validity region should be 
determined as follows. First of all, we have used assumption h>> w. Taking h = w as a criterium 
for the edge of our validity region one finds the dash-doted curves on figure |3l Strictly speaking 
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only those parts of solids curves are valid, which lie under the corresponding dash-doted curves. 
Now we have to ensure that rhombi chain is in quantum regime, i.e. currents are exponentially 
small. Simple estimates show that this is true above the doted line on figure 13 At large this 
limitation is very soft. Finally, our approximations are not valid at probabilities 'P^g very small 
or very close to unity since they rely on gaussian statistics for quantities A and B which certainly 
does not describe rare events. 

Qualitatively, results presented above can be interpreted as follows: at large Ej/Ec phase fluc- 
tuations in a single rhombi are weak and this makes random charges inefficient, however in a long 
chain global supercurrent is still exponentially suppressed. From the analysis presented above 
one concludes that we can guarantee with high probability existence of dominating, suppressed 
by quantum fluctuations, 4e-supercurrent in long chain {N ~ 20) with large Ej / Ec ^ 20 and for 
such a chain critical deflection 6^»'^/4>o ~ 10^^. If we compare this result with critical deflection 
(54>'^)^eg/^»o ~ 1.2 X 10^-^ for a regular chain with the same Ej /Ec^Q conclude that the influence 
of disorder on pairing effect is really small. On the other hand in a regular chain it was possible 
to obtain critical deflection of order 10^^ by choosing Ej/Ec ~ 8. For this set of parameters we 
cannot, use our theory quantitatively; we expect however that qualitatively the same behaviour as 
described above for larger values and EjjEc should be observed here as well. 



IV. MODULATION OF THE SUPERCURRENT BY CAPACITIVE GATE 

In this section we will discuss briefly influence of regular gates on pairing effect in maximally 
frustrated rhombi chain. We suppose that we have a gate, which is capacitively connected to all 
superconducting islands in one row (fig. [T]). 

Let us consider first rhombi chain with no disorder. Our system still can be described by 
equations (fT6lfT7b . but now a„ and hn are no more random. This coefficients can be expressed in 
terms of the gate voltage according to 

1 -|-cos7rC„y„ , 1 — cosTrCgV,, , 
an = ^ ^ ^ =a, bn = ^ ^ ^ = b (64) 

We will show here that the supercurrent in the chain is rather sensitive to the gate voltage. Since 
the gate voltage enters the problem in a particularly quantum way (via phases of virtual QPS), 
this dependence may provide an experimental test for the quantum nature of chain's state near the 
maximally frustrated point. 
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The problem we are considering now is much simpler as compared to the problem of random 
stray charges discussed above. Now our Hamiltonian commutes with the total spin of the chain 
and we can apply semiclassical approximation directly (e.g. by introducing spin coherent state 
path integral, c.f. It's easy to see that the ground state corresponds to the maximal total 

spin S = N/2. After proper redefinition of the time scale we can write imaginary-time action in 
the form 

5£ = Y J dx ^— zcos0(t) + -^^^ + (acos2xcos(|) + Z?sin2xsin(|)) sin 9^ = NSe{Ej /Eq, CgVg) 

(65) 

Here angles 6 and {|) parameterize coherent states of the spin. Again one easily find two types of 
tunneling trajectories corresponding to 4e- and 2e-supercurrents respectively (e.g. 4e-trajectory 
which connects |x = 0, = 7r/2, ^ = —n) with \x = n,Q = 7r/2, ^ = —n) and 2e-trajectory con- 
necting |x = 0, = 7r/2, (|) = —n) with \x = = 7r/2, ^ = 0)). Action Se is a function of 
Ej/Ec and CgVg and does not depend on the number of rhombi. 

From (l65t we have classical equations of motion (it is convenient to change variables according 
to0^7r/2 + /0) 

w 

— i:'-f (asin2jccos(|) — Z?cos2xsin(|)) cosh0 = (66) 
— (acos2xsin(|) — Z?sin2xcos(|)) = (67) 
cosh0(() — (acos2xcos(|)-|-Z7sin2xsin(|)) sinh0 = (68) 

Note that w/N ~ \) is a very good small parameter at Ej ::$> Ec- It allow us to obtain actions S^e 
and S2e analytically. Characteristic frequency of variable x is (Ox ^ ^/N/w ^ 1 , while for spins 
such a frequency is (O^ ~ 1. This means that on 4e-trajectory, which does not require any change 
in spin, we can assume spin to be constant (at least while coefficients a and b are not too close). 
This immediately leads to 



^4. ~ J—{\+cosnCgVg) (69) 



2^\ITj 

On the other hand, on 2e-trajectory smallness of w/N allows us to omit wx/N in equation for x. 
Resulting equations can be integrated analytically and we obtain 2e-action which, practically does 
not depend on Ej/Ec 



./|/j2_/,2| 2 J COSTlC gVg 

S2e ~ arcsh ^ = arcsh — — (70) 



1 —cosnCgVg 
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Figure 4: Modulation of 4e- and 2e-supercurrent in a clean chain by external gate. 



Actions and S^e can also be evaluated numerically. Results of numerical calculations are 
presented on figure IH Analytical results ( 16911701 are in very good agrement with numerical data. 
Note that 2e-supercurrent dominates over the current of pairs of Cooper pairs only in small vicinity 
of the point Cg V^/e = 1 . 

Note that in a regular chain change in 4e-action upon applying external gate voltage is nega- 
tive and proportional to the number of rhombi, i.e. external gate leads to significant increase of 
otherwise suppressed by fluctuations 4e-supercurrent. 

Let us now discuss influence of a regular gate on a chain with random stray charges. We will 
treat this problem within the same limit as in the previous section: we consider chain which is 
rather far from the maximally frustrated point in the sense h ^ w. From the preceding analysis we 
know that 4e-action is given by equation (l55t . Coefficient A acquires now the following form 



1 ^ 2 



(71) 



Here Ao is the value of A at zero external bias. Assuming that fluctuations of random charges are 
strong one easily finds 



(AA) = (AA^) = — ( sin^ TrCgVg + 8 sin' 



(72) 



So typical change in 4e-action can be estimated as 



55, 



Ae 



(73) 



This result differs from the analogous one for a regular (without offset charges) chain in two impor- 
tant aspects: i) we see that change in action is now proportional to y/N instead of A^^; ii) variation 
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of the tunnelling action can now be both negative or positive, i.e. applying external voltage we can 
decrease 4e-supercurrent as well as increase it, depending on the realization of charge disorder. 
From the above analysis we conclude that in a disordered chain current modulation by external 
gates is random and much weaker than in a regular chain, but it is still present and can be used to 
demonstrate quantum coherence of the tunnelling processes which occur in the chain. 



V. INFLUENCE OF "MAGNETIC DISORDER" UPON THE CROSSOVER POINT 

We now consider the effect of randomness in the values of magnetix fluxes through different 
rombi. Just as in the previous section we start from the grand partition function given by (l22b with 



dUn 



- (w COS 2x(x)5^ - hnS^^ Un (74) 

We presume that fluxes 4>" are uniformly distributed on (4>^ — A4>,4>r + A4>), i.e. probability 
distribution for hn is 

P{h) = <^ I U| ^^^^ 

I 0, \h-hQ\ > o 

where 

, 16A^(4>,-4>o/2) 16iVA4) 

ho = — , a = — — — . (76) 

Actually particular form of P{h) is not important for us as we assume that o <^ ho and use per- 
turbation theory in o/ho- We are interested in the critical deviation from the maximally frustrated 
point = 4>o/2 destroying 4e-supercurrent. Therefore we presume that ho ^ w. In such an ap- 
proximation we can use equation (l25l) to calculate Tri7„((3) and get effective action for variable x 

/■P x^(x) w^hN f 

S= dx^-^ — / cos2x(xi)D(xi-X2)cos2jc(x2) (77) 

Jo 2 A J 

1 ^ 

^(^i-'^2) = ;rT-rTi;exp(-/z„|xi-X2|) (78) 



At large N we can replace D(xi — X2) in dTTb by its mean value 

D(xi -12) = :^ [ dhe-^\^'-^^\p{h) 
2ho J 

Fourier transformation of D for small o/ho reads 



(79) 



S(a)) = _i_ ( i+a^ % ) (80) 
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Figure 5: Critical deviation from the maximally frustrated point in a chain with disordered flux in rhombi 
Again, using the fact that 

D-\(o) = oy' + hl + (5^- ^o'hl-j^ (81) 

we can perform Hubbard-Stratonovich transformation and find representation for partition func- 
tion with local action (after redefinition of time scale) 

S = /./rfx (f 4 4 + S) f + 1 +<vco.2.+ f-^y^ (82) 

Note that our formulation of the problem is a bit specific: we fix relative diversity o/Hq and look 
for ho which brings the chain to the point with equal Ae- and 2e-supercurrents. 

Discussion of the previous section concerning two types of tunneling trajectories and their 
connection to components of supercurrent can be literally applied to action (15^ . So we need to 
estimate classical action for 2e- and 4e-trajectories. Problem with a = was analysed in It 
was shown that 2e- supercurrent dominates at J ^ 1 while 4e- supercurrent — at J ^ 1. Crossover 
takes place at J ^ 3.24. From (l82b one can easily see that corrections to classical actions due to 
nonzero a should be of order (o/h)^. We analyse action (l82t numerically. Results are presented 
on figure 121 

It may come as a surprise, that with this formulation of the problem critical deviation 54>'^ 
(determined through mean flux in rhombi) grows with o/Ziq. However, this is quit reasonable 
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since zero deviation from the maximally frustrated point for one rhombi in a chain immediately 
prohibits 2e-supercurrent. Allowing for diversity of fluxes in rhombi we allow some of the rhombi 
to be closer to the maximally frustrated point than a "mean" rhombi. This fact causes strong 
suppression of 2e-supercurrent. 

So we conclude that randomness in fluxes is not important for the pairing effect, at least if the 
standard scatter of these fluxes A4> is not too large in comparison with the critical value S^)*^' found 
for a regular chain, cf-Eq- lfTTT) . 

VI. CONCLUSIONS 

In this paper we provide detailed calculations of superconductive current in a long frustrated 
rhombi chain with quenched disorder. We have considered two types of disorder: random stray 
charges and random fluxes in rhombi. We found that small (as compared to the critical deflection 
(54>'^)^g^ destroying 4e-supercurrent in a clean chain) fluctuations of fluxes piercing rhombi are 
not that important for the pairing effect. 

Main results of our paper concern effect of quenched random stray charges on pairing of Cooper 
pairs. For a chain which is relatively far from the maximally frustrated point in the sense ^ w we 
managed to calculate probability to find the system in the regime with dominating 4e-supercurrent. 
Stray charges, in principal, may significantly affect properties of the rhombi chain. In particular, 
in such a chain 2e-supercurrent exist even at the maximally frustrated point. However, as we have 
demonstrated in Sec. Ill, at large Ej/Ec and = 4>o/2 probability to find large 2e-supercurrent 
is small. This result itself is not a great surprise: in a perfectly classical chain stray charges have 
no effect at all. More important are two things: i) it is possible to combine low probability of 
finding significant 2e-supercurrent at the maximally frustrated point with exponential suppression 
of the supercurrent by quantum fluctuations; ii) in a perfectly classical chain critical deflection 
S^**^ scales with number of rhombi as 1/A^ (this can be easily seen from eq. © describing energy 
spectrum of classical rhombi chain) while in a disordered quantum chain this dependence is much 
weaker 54)^ ~ 1 /N^/^, c.f. 

In Sec. IV we have considered modulation of the supercurrent by capacitively coupled gate. In a 
regular chain, applying gate voltage one suppresses quantum fluctuations of rhombi and increases 
the supercurrent. Its dependence on the gate voltage is very strong: change of the logarithm 
of the supercurrent is linear in A^. On the contrary, in a chain with strong random stray charges 
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applying external gate can both increase or decrease supercurrent. Dependence of the supercurrent 
on the gate voltage is now much weaker: typical change of the logarithm of the supercurrent is 
now proportional to v^- Still we see that dependence of the current on gate voltage should be 
measurable. Such a dependence is one of the possible ways to demonstrate coherence of quantum 
phase slips in the chain. 
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